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1—^ , Abstract. I summarize recent results from lattice QCD on a few selected topics 

\f^ ' which are of interest to the heavy-ion physics community. Special emphasis is placed 

upon observables related to fluctuations of conserved charges and their connection to 

fH I event-by-event fluctuations in heavy-ion collision experiments. 
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Ph . 1. Transition temperature of QCD 



> 



At high temperatures and/or densities ordinary hadronic matter is expected to undergo 
1/^ I a transition to a new form of matter, viz. the Quark Gluon Plasma (QGP). The 

1^ I spontaneously broken chiral symmetry of Quantum Chromodynamics (QCD) is also 

O I expected to become restored at this point. The temperature at which such a transition 

t^^ I takes place is generally referred to as the QCD transition temperature (Tc). The value of 

Tc is one of the most basic input from lattice QCD to the heavy-ion phenomenology. At 
present almost all the precision state-of-the-art finite temperature/density lattice QCD 
computations are being performed using computationally cheaper staggered fermion 
/\ ' discretization scheme. At any non-zero lattice spacing (a) staggered fermions give rise 

c^ ' to a distorted hadron spectrum, which apart from containing a single Goldstone pion 

also contains 15 unphysical pions. Masses of these unphysical pions are larger than the 
Goldstone pion and the physical hadron spectrum is obtained only in the continuum 
limit a — )■ 0. Since the chiral property of the staggered fermions are sensitive to the 
masses of these unphysical pions it is essential to reduce their masses (relative to the 
mass of the Goldstone pion) while studying the QCD transition. 

One of the major progress made during the last couple of years in reducing this 
lattice artifact is the use of Highly Improved Staggered Quark (HISQ) action. The 
HISQ action largely reduces the masses of these unwanted pions even for relatively 
coarser lattices and consequently facilitates the continuum a — )■ extrapolation using 
moderate values lattice spacings. In this conference the HotQCD collaboration presented 
[1] their latest results on Tc obtained from observables related to the chiral symmetry 
restoration using two different types of fermion discretization {viz. HISQ and Asqtad 
actions) giving consistent results in the continuum a — )■ limit. HotQCD collaboration 
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Figure 1. (a) A sketch depicting the nature of QCD transition as functions of the hght 
(to„ — nid — rni) and strange (ms) quark masses at zero baryon chemical potential, 
(b) Evidence of 0{N) scaling of the order parameter as observed in Ref. [B]. See text 
for details. 



has also performed interpolations/extrapolations of their results to the physical values of 
quark masses using 0{N) scaling analysis. At the physical quark masses the continuum 
extrapolated value for the transition temperature is [H [2] — Tc = 157 ± 4 ± 3 ± 1 
(HotQCD preliminary). These latest results from the HotQCD collaboration is in 
excellent agreement with the previously published results (T^ = 147—157 MeV, obtained 
using observables which are sensitive to the light quarks) of the Wuppertal-Budapest 
collaboration [3l H] . 



2. Universal properties of the chiral transition 

The nature of the QCD transition depends crucially on the values of the quark masses. 



As for example, Fig. 1(a) shows a sketch of the nature of the transition as functions of the 



quark masses for a theory with two degenerate light (up and down) quarks with masses 
mi = iTLu = rrid and a heavier strange quark with mass rris at zero baryon chemical 
potential. In the limit m; — )■ and mg — )■ cxd, the relevant symmetry is isomorphic to 
the 3-d 0(4) spin model and the transition is second order belonging to the 3-d 0(4) 
universality class. For m^ = m^ — )■ the transition is first order. In the intermediate 
quark mass region a (rapid) crossover takes place from the hadronic to the QGP phase. 
The first order region is separated from the crossover region by a line of second order 
phase transitions belonging to the 3-d Ising Z{2) universality class. The first order 
region, the second order Z{2) and the second order 0(4) lines meet at a tri-critical 
point characterized by a certain value ml^ of the strange quark mass. While it is well 
established that for the physical values of the quark masses mi = mi- 
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the transition is a crossover [5] the location of the physical point in the context of Fig. 

^'^ then in the limit m^ — )• 



1(a) remains an open issue. More specifically, if m'^J^y'^ > m^' 

one should see a second order transition belonging to the 3-d 0(4) universality class. 

On the other hand, if m,^^^''^ < m}g then in m^ — )■ limit the transition should be first 

order. 

This question can be addressed by studying the the universal scaling behavior of the 



chiral transition as a function of decreasing light quark mass while keeping the strange 
quark mass fixed at its physical value. In the vicinity of a second order phase transition 
the behavior of the free energy of a system is largely governed by the universal scaling 
properties— / {mi, rris, T, (xb) = h^^^/^fs{z)+freg {mi, m^, T, (ib) , where h = mi/{homs) 
and t = {T — To)/{toTo) are, respectively, the explicit (chiral) symmetry breaking field 
and the reduced temperature variable. Tq denotes the critical temperature in the chiral 
limit {h — > 0) and z = th'^^^^ is the so-called scaling variable. The scaling properties of 
the order parameter (~ df /dh) is given by M = nis (^\I') /T'* = h^/^ fG{z), where {|^\l/) 
is the chiral condensate. The critical exponents (/3, 5) and the scaling functions {fs{z), 
fG{z)) uniquely characterize the universality class of the transition. The parameters 
to, hQ and To are non- universal, i.e. they depend on the action, lattice spacing, value 
of the bare strange quark mass etc. . The function freg{mi,ms,T, (ib) is not related to 
universal scaling properties and depicts the usual non-critical regular part of the free 
energy. 

Recently such scaling studies were performed in Ref. [6] by the BNL-Bielefeld 
collaboration using p4fat3 staggered quark action. These studies show that in the limit 
rrii — ;■ 0, keeping rng = mPj^'^'^, the chiral transition belongs to 3-d 0{N) universality 
class Ij. As for example. Fig. 1(b) shows that the properly scaled order parameter. 



h~^/^M, agrees with the relevant 0{N) scaling function /g(^) for nii < 771^/20 which 
corresponds to a Goldstone pion mass of m^ ~ 150 MeV. These studies indicate that 
the m^^^''^ is probably larger than ml^. Furthermore, they also indicate that the physical 
QCD {rrii ~ ma/27) may also lie within the scaling regime of the actual chiral transition. 
However, these studies where performed using p4fat3 staggered action which is known to 
have large discretization errors. Hence, in order to confirm these results similar studies 
have to be performed using the HISQ action. Such studies are currently underway. 

3. Curvature of the phase transition line for small /i^ 

Based on the indications from the scaling studies mentioned in the previous section 
the conjectured phase diagram of QCD in the limit of massless light quarks is depicted 



in Fig. 2(a) At zero chemical potential (/i^) the chiral transition is second order, 
belonging to the 3-d 0(4) universality class. This 0(4) chiral phase transition line 
extends in the T — fiB plane and meets the 1st order phase transition line emanating 
from the T = axis at the so-called QCD (tri-) critical point which belongs to 3-d Z{2) 
universality class. The universal scaling properties of the chiral transition has been 
extended to compute the curvature of the 0(4) chiral critical line for small values 
of baryon chemical potential Hb/T « 1. Since the baryon chemical potential does 
not explicitly break the chiral symmetry, from the perspective of the universal scaling 
properties of the chiral transition, /i^ is equivalent to the thermal variable. Further, 

X As mentioned before in Section [1] for staggered fermions at non-zero lattice spacings there is only one 
Goldstone pion in the chiral limit as opposed to three in the continuum. Hence, for staggered fermions 
the relevant symmetry group is 3-d 0(2) instead of the 3-d 0(4) of the continuum. 
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Figure 2. (a) A sketch of the conjectured phase diagram of QCD in the temperature- 
chemical potential plane in the limit of massless light quarks, (b) Scaling analysis of 
the "thermal" susceptibility of the order parameter and determination of the curvature 
of the chiral critical line [^ . (c) Comparison of the phase boundary in the T — \ib 
plane and the chemical freeze-out curve. See text for details. 



due to CP-symmetry of the free energy, the reduced temperature variable must be even 
in [iB — ^ = [(^ — ^o)/^o + ^B {,\^bIT^ ] Ao- Here we assume [IbJT << 1 and consider 
only the leading order term in ^b/T. The 0(4) chiral transition line inT — ^b plane is 
defined by the condition t = and hence at the leading order in fiB the chiral transition 
temperature is given by — To(/iB)/To(0) = 1 — KB[fJ'B/To{0)]'^. Thus kb is the curvature 
of the chiral transition line around fiB = 0. The value of the this curvature parameter kb 
can be determined by studying the universal scaling properties of the so-called "thermal" 
susceptibility of the order parameter — Xm,B = d'^M/d{^,B/TY = 2K,Bh^'^~^^^^^ fQ{z)/to, 
where /^(-z) is a scaling function uniquely associated with the specific universality class. 
Once the non-universal scaling parameters to, Hq and Tq{0) is determined from the 
scaling analysis at /is = 0, as have been done in Ref. [6], scaling analysis of Xm,B can 
be preformed in order to determine the only unknown parameter Kb- 

Such a scaling analysis of the "thermal" susceptibility of the order parameter was 
performed by the BNL-Bielefeld-GSI collaboration in Ref. [7], using p4fat3 staggered 
fermions with two different lattice spacings, and the curvature of the chiral transition line 
was determined to be kb = 0.0066(9). In Fig. 2(b) we show the scaling analysis of the 
"thermal" susceptibility of the order parameter from Ref. |§|. Once the the curvature 
of the chiral transition line is known it is possible to compare the phase boundary in 
T — fiB plane, given by To(/iB)/To(0) = 1 — HBlfJ'B /To{0)]'^ , with the (chemical) freeze-out 
curve [8j. This gives an idea about how far the chiral critical line is from the freeze-out 
of the heavy-ion collision experiments. Such comparison is shown in Fig. 2(c) In the 



limit rri; —7- oo, i.e. for a pure gluonic theory, the transition temperature is independent 
of hb and in the limit mi ^ the phase boundary is the one obtained by the BNL- 
Bielefeld-GSI collaboration in Ref. [7]. The transition line for the physical quark mass 

§ In Ref. [7] the "thermal" susceptibility of the order parameter Xm.,q was defined with respect to the 
average light quark chemical potential /ig = {^u + Md)/2. Since for a two flavor theory /i^ = M-b/3, 
the curvature parameter Kg of Ref. [7] is trivially related to the curvature parameter discussed here 
through the relation kb = Hq/9. 
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Figure 3. (a) Second order susceptibility for the baryon number, (b) fourth order 
susceptibility for up quark number and (c) sixth order susceptibility for the electric 
charge. The ne"w calculations for the HISQ action have been performed by the BNL- 
Bielefeld collaboration and that for the stout action are from the Wuppertal-Budapest 
collaboration [TT]. The older data for the p4fat3 action is taken from Ref. [T^]. The 
dotted lines indicate the values of these quantities for the hadron resonance gas (HRG) 
model and for the free fermionic gas (SB). 



should lie in between these two limits. In fact, the phase boundary for the physical quark 
masses has also been computed recently from lattice QCD by the Wuppertal-Budapest 



collaboration [3l |9]. As can be seen from Fig. 2(c) the phase boundary for the physical 
quark mass almost coincides with that in the chiral limit and they agrees with the 
freeze-out curve [8] till hb ~ 150 MeV and then slowly starts to deviate. However, one 
has to remember that the lattice computations for the phase boundary are at present 
performed only upto the leading order in /ig. 



4. Fluctuation, correlations: lattice QCD meets experiments 

Although lattice QCD has proven to be a successful technique for performing ab- 
initio, quantitative studies of QCD a direct computation at non-zero baryon chemical 
potential is, unfortunately, plagued by the infamous sign problem. However, important 
information at non-zero chemical potentials may be obtained using the well established 
numerical technique of Taylor expansion of the QCD partition function with respect 
to the baryon chemical potential [10]. In this method one Taylor expands the 
partition function (Z)/free energy (/) /pressure (P) in power series of the /xs 
around /i^ = 0- P{^ib,T)/T^ = Er=o(V^!)x^(T)(/iB/T)", where VT\^{T) = 
[d"'Z{fiB,T)/d{^B/T)"']^g=o and V is the volume of the system. Xn ^^^ ^^e generalized 
baryon number susceptibilities and since they are defined at /ib = usual lattice QCD 
simulations at zero chemical potential can be used to compute these susceptibilities. 
Similar susceptibilities Xn^ Xn ^'tc. can be defined from the Taylor expansions of the 
pressure with respect to the electric charge chemical potential fiq , the strangeness 
chemical potential fis etc. . 

Since the generalized susceptibilities are derivatives of the free energy with respect 
to the chemical potentials universal scaling properties of these susceptibilities can 
also be readily predicted. From discussions in Sections [2] and |3] it can be inferred 
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that two derivatives of the free energy with respect to fiB is equivalent to a single 
derivative with respect T and consequently the scaling behavior of the susceptibilities 
are — Xn ~ /^(i-a-"./2)yw >(^2:) + Cn,reg- Here /("/^)(z) are unique scaling functions 
and Cn,reg ^re contributions from the non-critical regular part of the free-energy. In 
order to understand the universal scaling behavior of the generalized susceptibilities 
it is also important to note that for 3-d 0(4) universality class the critical exponent 



— 1 < a < 0. In Fig. 3(a) we show the second order baryon number susceptibility. Close 
to Tc, X2 ~ C2,reg T A± \{T — Tc) /Tc\ ~ is uou-divergeut and behaves like the energy 
density. While for T < T^. the new lattice data is in good agreement with the hadron 
resonance gas model, for T > T^ approach towards the ideal gas value is much smoother 
and slower compared to the older data from the p4fat3 action [12]. Also note the good 
agreement between the new data coming from the HISQ (BNL-Bielefeld collaboration) 
and the stout [H] actions. Close to Tc the behaviors of the higher order susceptibilities 
are also consistent with there expected scaling behaviors. The fourth order susceptibility 



(see Fig. 3(b) ) X4 ~ C^^reg ± ^± | (T — Tc)/Tc\~°' remains non-divergent and behaves like 
the specific heat by always staying positive and showing a cusp like structure in the 
vicinity of Tc. The sixth order susceptibility xe ~ 1^A±\{T — Tc)/Tc\'^ "^ changes 



sign across T^. (see Fig. 3(c)) and diverges at T^ in the chiral limit. Based on the later 
observation that even for /i^ ?^ Xe < in the vicinity of T^, it has been recently argued 
[13] that if the measured values of the sixth order moments of the conserved charges 
come out to be negative for the on-going heavy-ion collision experiments at the Large 
Hadron Collider (LHC) then it may indicate that the transition from the hadronic to 
the QGP phase happens very close to the freeze-out curve. 

The off-diagonal susceptibilities xfi^; Xii ^'^c. , are measures of correlations among 
different conserved charges. As a result these quantities are good probes for the degrees 
of freedom inside QGP. As for example, inside QGP where strange quarks {B = 1/3, 
S = —1) are the relevant degrees of freedom strangeness can only exists in direct 
conjunction with the baryon number. On the other hand, at very low temperatures in the 
hadronic phase where the relevant strangeness carrying excitation are the kaons [B = 0) 
there is no apparent correlation among the baryon number and strangeness. Based on 
these observations the observable baryon strangeness correlation Cbs = — 3xfi^/X2 ^^^ 



proposed in Ref. [11]. In Fig. 4(a) we show the this quantity computed using the 
HISQ (by BNL-Bielefeld collaboration) as well as using the stout [11] action. Below T^ 
this quantity is in very good agreement with the results predicted by the HRG model 
indicating that at these temperatures the strangeness is truly carried by the hadrons. 
On the other hand, for T > 250 MeV the baryon strangeness correlation is quite close to 
its ideal gas value indicating that above these temperatures the strangeness and baryon 
number are carried by the strange quarks. Also note that the lattice data from both 
the actions are in excellent agreement with each other. 

So far we discussed the physics of the generalized susceptibilities at zero chemical 
potentials. However, susceptibilities for non-zero values of chemical potentials have more 
direct relevance for the heavy-ion collision experiments. These can be obtained from 
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Figure 4. (a) The baryon-strangeness correlation [T3] using the HISQ (BNL-Bielefeld) 
and the stout [11 actions. The dotted lines indicate the results for the hadron 
resonance gas model (HRG) and for the free fermionic gas (SB), (b) Experimentally 
measured moments of net-proton fluctuations |15) are compared to lattice QCD 
computations of ratios of baryon number susceptibilities by the BNL-Bielefeld [IT] 
and the Mumbai [TB] groups. The HRG results are from Ref. [TO], (c) An example of 
the sensitivity of the ratios of conserved charge susceptibilities to different orders of 
Taylor expansion [17]. See text for detail. 



the Taylor series expansions of the susceptibihties themselves, e.g. the baryon number 
susceptibilities at /i^ > are given by— Xnil^-B^T) = Efclo(l/^0x^+fe(O,T)(/iB/T)'=. 
These generalized susceptibilities are measures of the moments of the conserved charge 
distributions. As for example— VT^xf{f^B,T) = (Nb), l/T^xf (/iB,T) = {{6NBf), 
VT\^{liB,T) = {{5NBf), VT\f{^iB,T) = {{5NbY) -?,{{5NBff etc. measure 
different moments of the distribution of the net baryon number NBifJ^BjT), where 
5Nb{^ibiT) = Nb{iJ'B,T) — {NB{ftB,T)). On the other hand, such moments of 
conserved charges have also been measured by the STAR experiment [15] using event- 
by-event fluctuations in the Relativistic Heavy Ion Collider (RHIC) at various energies 
{^/SNN) ■ As for example, in heavy-ion collision experiments one measures the 
moments — mean MB{y/sNN) = {Nb), variance (J%{^Jsnn) = {(SNb)'^), skewness 
SB{^/sNN) = {{^NbY) /a%, kurtosis Kb{^snn) = {{^NbY) /a% - 3 etc. . Assuming— 
(i) experiments measure conserved charge distributions of a thermalized system, (ii) 
measured moments characterize the chemical freeze-out condition and (iii) T{y/SNN) 
and fiB{y/S]\fiy) at the chemical freeze-out can be modeled by the hadron resonance 
gas model |8j, the experimentally measured volume independent combinations of 
the moments can be related to the ratios of susceptibilities computed from lattice 

QCD— (Tl{y/SNN)/MB{^ySNN) = X2 if^B,T) / xf (fJ^ByT) , CTfi ( a/s^) Sfi ( ^S^^) = 
xiil^B,T)/x2ifJ'B,T), al{y/SNN)KBiy/SNN) = xf il^B,T) /X2 ifJ-B^T) ctC. . 

As net proton fluctuations can be treated as a proxy for the net baryon number 
fluctuations [I6], in Fig. 4(b) we show a comparison of the experimentally measured 



moments of net proton distribution [15] to the lattice QCD calculations of the ratios of 
the baryon number susceptibilities using improved p4fat3 [IT] and un-improved naive 
[18] staggered actions. As can been seen, the experimental data and the lattice QCD 
computations are in good agreement with each other and also with that from the hadron 
resonance gas calculations [T^]. This, for the first time, shows a direct comparison 
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between the heavy-ion colhsion experiments and lattice QCD computations. However, 
in order to make such comparisons it is essential to perform reliable computations of the 



higher order susceptibilities. As for example, in Fig. 4(c) we show the results [17] for 
ratio xf /xf computed in leading and next to leading orders of the Taylor expansion. 
As can be seen, around ^/SNN ~ 20 GeV (i.e. fiB ~ 200 MeV) the contribution of the 
next to leading order term is much larger than the contribution of the leading order 
term. Moreover, at present these lattice results are not continuum extrapolated and 
were calculated using actions which are known to suffer from significant discretization 
effects and resulting in larger values for the transition temperature. Hence, in future, 
it is essential to carry out these lattice computations using HISQ like highly improved 
actions with accurate determination of the higher order susceptibilities. This is specially 
important for the observables like moments of net electric charge fluctuations [15] as 
they are sensitive to light charged hadrons like the pions (see discussion in Section [1]). 
Such calculations are currently in progress 
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